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An elementary proof is given of the fact that if A, B are nonsingular nX n matrices over a principal 
ideal ring /?, then the kth invariant factor of AB is divisible by the kth invariant factor of A and by the 
kth invariant factor of B, 1 ^ A ^ n. Some consequences are drawn. 
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1. Introduction 

Let R be a principal ideal ring (a commutative ring with identity 1 in which every ideal is 
principal). If AeR n (the ring of n X n matrices over R) A T will denote its transpose. If in addition 
A is nonsingular then 

S(A) = diag( Sl (A),s 2 (A), . . . 9 sn(A)) 

will denote the Smith normal form of A (see [2] for an excellent reference on this topic). It is well- 
known that if A, B are nonsingular elements of R n then the determinantal divisors of AB are divisible 
by the corresponding determinantal divisors of A and of B. It is not so well-known that the same 
result is true for the invariant factors: i.e., Sk(AB) is divisible by Sk(A) and by Sk(B), 1 ^ k ^ n. 
An interesting consequence is that S(AB) =S(A)S(B) , provided that A, B have relatively prime 
determinants. This result is a consequence of a rather general theorem about rings which is given 
by Kaplansky in his paper [l]. 1 Since Kaplansky did not include a proof of his theorem in his paper, 
and since the proof at any rate would be ring-theoretic, it is desirable to have a purely elementary 
matrix-theoretic proof, and this is what is given here. 

We also note that the result concerning the determinantal divisors follows as a corollary, 
since if A is nonsingular then its A;th determinantal divisor dk(A) is just s 1 (A)s2(A) . . . Sk(A) , 
l^k^n. 

2. A Lemma 

We first prove the following lemma: 

Lemma 1: Suppose that I I is a nonsingular element ofK n which is in Smith normal form. Let 

m be any non-zero element ofK and suppose that there is a matrix U ofK n such that (det U, m)= 1, 
and 

(TT A \ /* *\ 

0KrW modm - 
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Then K = mod m. 
Proof: Put 

H = Amg (hi, h 2 , . . ., h r ), K—diag (k\,k 2 , . . -,k s ), 

where r + s = n and hi\h\+\ (l^i^r— 1), kj\kj+\(\ ^y^s — 1), h r \k\. Partition U so that 



^r 1 



v 3 c/ 4 ;' 



where f/i is r X r, f/ 4 5 X 5. Then (1) implies that 

(2) U 3 H = mod m, 

(3) U 4 K=0 mod ra. 

We can multiply (2) on the right by diag (h r /h u h r \h 2 , . . ., 1) to obtain 

/&/•£/:$ = mod m. 

Put 

(4) (fc r , m) = 8. 



Then 



and since (A r /8, ra/8) = 1, 



o o 
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(5) f/ 3 =0 mod - 

o 

Now set K=k\K\ where #' = diag (1, k 2 \k\, . . ., k s /k\). Then from (3), 

kiU 4 K f = mod /n. 
Put 

(6) (*,,ro)=A. 
Then as before we deduce that 



(7) UJC^O mod^- 

Now 8\ A, in virtue of (4), (6), and the fact that h r \ki. It follows that m/A| m/8, and so (5) holds modulo 
ml A as well. Thus 

TJ (U x U 2 \ (Ut U t \ A m 

U= \U 3 Uj = \0 Uj m ° d A^ 

det U = det U t det Ut mod -r • 
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Since (det £/, m) =1, it follows that (det C/ 4 , ra/A) = 1. Hence (7) implies that /T =0 mod m/A, 
and so A = mod ra, since the 1,1 element of K' is 1. Thus (6) implies that m\k\ , and the conclusion 
follows. 

3. The Theorems 

We are now prepared to use Lemma 1. Let A, B be nonsingular elements of R„. Then matrices 
U, C of R n exist such that U is a unit matrix and 

(8) US(AB)=S(A)C. 

Certainly s\(A) divides S\(AB), since S\(A) is the greatest common divisor of the elements of 
A and s\ (AB) the greatest common divisor of the elements ofAB. For 2 ^ k ^ n— 1, choose m = 5a (/4 ) 
and apply Lemma 1. We are left with k = n. Write E/= (^/j), C= (cy). Then (8) implies that 

u ij Sj(AB)=c ij Si(A), 

so that for i = n, 

u tl jSj(AB) =0mod5 /( (^). 

Since Sj(AB) \s„(AB) , l^j^n, this implies that 

u n jS n (AB) = mod s,, (.4 ) . 

The fact that 5/, (A) \ s,,(AB) now follows, since (u u \, u, r >, . . .,u wn ) = l. 

If we note that in addition S(A T )=S(A), the entire argument may also be applied to the 
pair B l \A l \ and we finally obtain 

THEOREM 1: Let A, B be nonsingular elements o/R n . Then s k (AB) is divisible 6ys k (A) and by s k (B) 
/or 1 ^ k ^ n. 

From this theorem we easily deduce 
THEOREM 2: Suppose that A, B arc elements o/R n with relatively prime determinants. Then 

S(AB) = S(A)S(B). 

PROOF: Since (det A, det B) = l and s k (A) | det A, s k {B) | det #, it follows that (s k {A), s, : (B)) = l 
for 1 ^ A ^ n. Then Theorem 1 implies that 

s k (AB) =0 mod s k (A)s k (B) 9 1 ^ A: ^ /i. 

// n 

But J| Sk-(AB) is a unit multiple of det(/4#), and J| s k (A)s k (B) is a unit multiple of det A det #. 

A = 1 fc= 1 

It follows that s k {AB)\s k (A)s k (B) is a unit for 1 ^ k ^ n. But this implies that in fact s k (AB) = 
s k {A)s k (B) for 1 ^ k ^ n, since associated elements in corresponding diagonal positions of matrices 
in Smith normal form must be equal. This completes the proof. 

4. Concluding Remarks 

Theorem 2 is definitely false if (det A, det B) > 1. Thus if /tt is any element of R, and 

h; :.)• Hi. :> 

then 
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S(AB) 



but 

'm 
\0 m J' 

This example also shows that S(AB) need not equal S(BA), since here 

1 



However, S(AB) and S(BA) are equal if ^4 and B have relatively prime determinants as is readily 
seen from Theorem 2. 

A simple example illustrating the use to which Theorem 1 may be put is furnished by choosing 
A as the incidence matrix of a finite projective plane of order n, so that A is a v X v rational integral 
matrix satisfying 

AA T = A T A = nI + J, 

where v = n 2 + n + 1 and J is the matrix all of whose entries are 1. It is easy to show that the 
invariant factors of nl -h J are 

l(once), n(v — 2 times), n(n + l) 2 (once). 

Thus if the Smith normal form of A is 

S(A) = diag (ai,a 2 , . • . , «v), 
then 

«x= 1, «i|n(2 ^ i *s= V- 1), a^l/iC/i-f l) 2 . 

Now the facts that 

v-l 

aia 2 . . . a v = n 2 (n-\- 1) , 
and (a*, n + 1) = 1 for 1 ^i^?; — 1 , imply that a y = (n + l)ai, where now 



aia 2 . . . aiy-ia^ = n L , 




ai|ai+i(l ^ i ^ w — 2),a»_i|ai,ai|n. 




Choosing ra square-free, we easily obtain 




COROLLARY 1: Let A be the incidence matrix of a finite projective plane of order n, where n is 




square-free. Then the invariant factors of A are 




/n 2J rn . \ ( n 2 + n . \ w x 
11 -hi tames 1, /z 1 - 1 tawies 1, ai(ai -r 1) (orace). 




Of course such a matrix is known to exist only if n is a prime, and this result might possibly 




be of some use in settling the question of existence for other square-free values of n. 
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